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ABSTRACT 

We impose constraints on the topology of the Universe determined from a search for 
matched circles in the temperature anisotropy patterns of the 7- year WMAP data. We 
pay special attention to the sensitivity of the method to residual foreground contami- 
nation of the sky maps, and show that for a full sky estimate of the CMB signal (the 
ILC map) such residuals introduce a non-negligible effect on the statistics of matched 
circles. In order to reduce this effect, we perform the analysis on maps for which the 
most contaminated regions have been removed. A search for pairs of matched back-to- 
back circles in the higher resolution WMAP W-band map allows tighter constraints 
to be imposed on topology. Our results rule out universes with topologies that predict 
pairs of such circles with radii larger than Ofmin ~ 10°. This places a lower bound on 
the size of the fundamental domain for a flat universe of about 27.9 Gpc. This bound 
is close to the upper limit on the size of Universe possible to detect by the method of 
matched circles, i.e. the diameter of the observable Universe is 28.3 Gpc. 
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1 INTRODUCTION 

According to General Relativity, the pseudo-Riemannian 
manifold with signature (3,1) is a mathematical model of 
spacetime. The local properties of spacetime geometry are 
f"*) described by the Einstein gravitational field equations. How- 
^-H . ever, they do not specify the global spatial geometry of the 
IL" ' universe, i.e. its topology. This can only be constrained by 
. 5^ ' observations. 

I The concordance cosmological model assumes that the 
5^ ' universe possesses a simply-connected topology, yet recently 
detected anomalies on large angular scales in the cosmic 
microwave background (CMB) anisotropy suggest that it 
may be multiply-connected. Evidence of such anomalies 
comes from the suppression of the quadrupole moment, 
alignment of the quadrupole and octopole and asymme- 
try in the temper ature anisotropy observed in two hemi- 
spheres on the sky ( de Oliveira-Costa et al.|[2cio3; Copi et al.l 



iKunz et aO |2006|. I2OO8I). and t h e distribution of matched 
patterns (ICornish et al.l |2004|: Aurich. Lustig. fc SteineJ 



l2005l . l2006inThenll2006l : iKev et al.ll2007l ) 



20041 : lEriksen et al 
20o4 l 



20041 : [Hansen et al]|2004l : ISchwarz et"al] 



The present constraints on topology were placed by 
studying two signatures of multi-connectedness in the CMB 
maps: the large scale damping of the power in the direction 
of a small dimension of the domain, which causes a break- 
down of statistical isotropy l|de Oliveira-Costa et al.l |2004| : 



In this work, we constrain the topology of the Uni- 
verse using the method of matched circles proposed by 
Cornish. Spergel. fc Starkman (il998 j l a nd apply it to the 7- 
year WMAP data ( Jarosik et al. 201Cf ). In contrast to the 
majority of previous studies, we will pay special attention 
to the impact of Galactic foreground residuals on the con- 
strain t s. So me consideration of this problem was made by 
iThenj (|2006l ') in their analysis of the first year WMAP data 
release. We will use also Monte Carlo (MC) simulations for 
the estimation of the false detection level of the statistic. 
The method is applied to higher resolution maps than pre- 
viously, which implies a lower level of false detection and 
therefore tighter constraints on the size of the Universe. As 
a result of computational limitations, we will restrict the 
analysis to a search for back-to-back circle pairfl. 

In the following two sections we describe data used in 
analysis and simulations of the CMB maps for a flat uni- 
verse with the topology of a 3-torus which were used to test 
the reliability of our implementation. The statistic adopted 
in our studies is presented in Section |4] The results and 
conclusions are described in the last two sections. 
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2 DATA 

The search fo r matched circles wa s performed on the 7-year 
WMAP data Jjarosik et al.ll20ld ). The WMAP satelhte ob- 
serves the sky in five frequency bands denoted K, Ka, Q, 
V and W, centred on the frequencies of 22.8, 33.0, 40.7, 
60.8, 93.5 GHz with angular resolutions of 52.8', 39.6', 30.6', 
21' and 13.2', res pectively. The map fl are pixelized in the 
HEALPI5<0 scheme (|G6rski et al.ll2005l ) with a resolution pa- 
rameter A'^side = 512, corresponding to 3 145 728 pixels with 
a pixel size of ~ 7 arcmin. 

In particular we used the Internal Linear Combination 
(ILC) map and the W-band map corrected for Galactic emis- 
sion outsid e of the applie d mask using a template fitting 
scheme ((Gold et al.l I2OI0I ) and smoothed with a gaussian 
beam profile of full width at half-maximum (FWHM) — 20' 
to decrease noise level. The former has been used for a num- 
ber of full sky analyses, and versions available with earlier 
data releases have already been exploited in the search for 
matched circles. The latter is the highest angular resolu- 
tion data measured by the WMAP satellite. Since the false 
detection level for such a map is lower, it provides the op- 
portunity to estimate tighter constraints on topology. We 
electe d not to use data preprocessed as in ICornish et al.l 
l|2004h . i.e. a noise-weighted combination of the Q, V and 
W-band maps smoothed to Q-band map resolution outside 
of the Kp2 Galactic cut, and the ILC map within this re- 
gion. Such a map has complex properties since it consists 
of two regions with different angular resolutions; moreover, 
as we will see the ILC map remains too contaminated by 
Galactic foreground residuals inside the Kp2 cut to be used 
for cosmological analysis. 



3 SIMULATIONS OF THE CMB 
ANISOTROPY MAPS IN THE 
MULTI-CONNECTED UNIVERSE 

To test the reliability of codes used to search for the sig- 
nature of a multi-connected topology in the Universe, we 
performed simulations of CM B skies for a flat un iverse 
with the topology of a 3-torus (|Riazuelo et al.ir2004al V The 
simulations can be made in two ways: one can compute 
the covariance matrijfl of the sphe rical harmonic coeffi- 
cients aim, M. ern.(em)' = {o«m'J£'m')i (|Riazuelo et al.|[2004bl : 
iPhillips fc Kog ut 2006) and then generate correlated co- 
efficients using the Choleski decomposition of the covariance 
matrix M = L L^ and uncorrelated unit variance Gaussian 
variable a;, a = L a;, or one can directly compute the aim 
coefficients using equation 



k \k\=k 

(1) 

where Ar = tq — Tr is difference between the present con- 
formal time To and the time of recombination Tr, P{k) is 
the density perturbation power spectrum, (fc, At) is the 
response function for temperature, Ek is a Gaussian random 

^ available at http://Iambda .gsfc. nasa.gov| 

^ http://healpix.jpl.nasa.gov 

^ assuming zero mean of the CMB maps 



variable which satisfies (efcej/) — p^yj-^fefe' and V is the vol- 
ume of the fundamental domain. The sum over wavevectors 
k was split into two sums "^^^ = X^fc X]|te|=fc " f^e sum over 
all possible magnitudes of the k vector and over all fc of 
magnitude fc. 

The advantage of the latter approach is a less aggres- 
sive scaling in the computation time, ^nmx't^ (number of aim 
coefficients oc ^max times number of wave modes oc ^^ax^) 
in comparison with the formei[f| that goes as i^^xV- On the 
other hand, the disadvantage of the latter approach is that 
for each simulation we have to repeat the computation^. 
Nevertheless, we decided to adopt this method since the 
difference in scaling is particularly significant for the high 
resolution maps of interest in our studies. We also do not 
need more than one map for a given topology to evaluate 
the possibility of detection. 

The topology does not affect local physics, so the equa- 
tions governing the evolution of cosmological perturbations 
are left unchanged. Thus, to compute the response func- 
tions, A£^{fej_AT^}j_one_£an_jise_tlie^ available CAMB 
code (|Lewis. Challinor. fc Lasenbvl200ol ) . A change of topol- 
ogy translates into a change of the modes that can exist in 
the universe. Therefore, we need to compute the response 
function only for a set of allowed wave numbers fc. To speed 
up computations, we use in our code precomputed response 
functions. Then, the most time consuming part is the com- 
putation of the sum X]|k|=fc Yl^{flk)ek. 

The simulations do not include gravitational lensin g of 
the C MB as the lensing deflection angle is small ( Seli ak I 
Il996l ). ~ 3', compared to the minimal angular scale taken 
into account in the simulations, ~ 20'. The effect of the fl- 
nite thickness of the last scattering surface is included in the 
simulations through the response function, A|^(fc, Ar), com- 
puted by CAMB. However, the effect starts to be important 
at angular scales smaller than ^ 4' and will therefore not be 
studied here. 



3.1 Requirements and numerical implementation 

To study the signatures of a given topology, a CMB map 
is required with resolution comparable to the angular size 
of the beam profile used to smooth the WMAP W-band 
data (i.e. ~ 20') to be analysed here (see Section [2}. We 
have adopted ^max = 500 in our simulations. The dimension 
of the fundamental domain of the 3-torus was L = 2c/ Ho, 
which is about three times less than diameter of the observ- 
able Universe i.e. ~ 6.6c/ Ho- In such a universe there are 
many pairs of matched circles of different radii. The time 
needed for the generation of one such CMB map on a single 
processor with clock speed 3 GHz is about 42 hours. 



^ the parity and symmetry relations can reduce by an order of 
magnitude the number of different elements in the matrix and 
computation time. 

^ in the former case one need only generate a vector of Gaus- 
sian random numbers and multiply by the Choleski decomposed 
matrix. It is much faster than the full computation. 
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4 SEARCHING FOR THE CIRCLES IN THE 
SKY 

If light had sufficient time to cross the fundamental cell, 
an observer would see multiple copies of a single astronom- 
ical object. To have the best chance of seeing 'around the 
universe' we should look for multiple images of the furthest 
reaches of the universe. Searching for multiple images of the 
last scattering surface - the edge of the visible universe - 
is then a powerful way to constrain topology. Because the 
surface of last scattering is a sphere centred on the observer, 
each copy of the observer will come with a copy of the last 
scattering surface, and if the copies are separated by a dis- 
tance less than the diameter of the last scattering surface, 
then they will intersect along circles. These are visible by 
both copies of the observer, but from opposite sides. The 
two copies are really one observer so if space is sufficiently 
small, the CMB radiation from the last scattering surface 
will demonstrate a pattern of hot and cold spots that match 
around the circles. 

The idea of using these circle s to st udy topology is due 
to ICornish. Spergel. &: StarkmanI (|l998l ). Therein, a statis- 
tical tool was developed to detect correlated circles in all 
sky maps of the CMB anisotropy - the circle comparison 
statistic 



Sp,r{0t, 



= 



(2) 



where {) = J^^ drf) and Tp{±(j)), Tr{(j) + (j>,) are the tem- 
perature fluctuations around two circles of angular radius a 
centered at different points, p and r, on the sky with relative 
phase <j>,. The sign ± depends on whether the points along 
both circles are ordered in a clockwise direction (phased, sign 
-|-) or alternately whether along one of the circles the points 
are ordered in an anti-clockwise direction (anti-phased, sign 
— ). This allows the detection of both orientable and non- 
orientable topologies. For orientable topologies the matched 
circles have anti-phased correlations while for non-orientable 
topologies they have a mixture of anti-phased and phased 
correlations. The statistic has a range over the interval 
[—1, 1]. Circles that are perfectly matched have S = 1, while 
uncorrelated circles will have a mean value of S — 0. To 
find matched circles for each radius a, the maximum value 
'S'max(Q^) ~ max S'j^^(Qf, (^*) is determined. 

In order to speed up the computations, one can 
use the fast Fourier transform (FFT) along the cir- 
cles, Tp[(j>) = 5]]^ Tp^m exp(im(/!>), and, by rewriting the 
statistic as Sp^r{o-,4>*) = exp(— imi;/)*), where 

8m = 2X:^Tp,^T*™/E„(|Tp.n|' + |rr,„n, use iuversB 
fast Fourier transform. It reduces the computational cost 
SjJ;r(a, (f)*) from ~ A'' to ~ A''"'^''^ log A*' operations, where 
is number of pixels in the map. 

In our studies w e will use the statistic as modified in 
ICornish et al.l (|2004l ) to include weighting of the mth har- 
monic of the temperature around the pth circle Tp,m by the 
factor |m| that takes into account the number of degrees of 
freedom per mode: 



Sp,r{a,<i>*) 



EJ^I(|rp,„p + |T.,„p) 



(3) 



where the Fourier coefficients Tp,™ are complex conjugated 
in the case of the Sp^^ statistic. Such weighting enhances the 



contribution of small scale structure relative to the large 
scales fluctuations. It is especially important because the 
large scale fluctuations which are dominated by the inte- 
grated Sachs- Wolfe (ISW) effect. This can obscure the im- 
age of the last scattering surface and reduce the ability to 
recognise possible matched patterns on it. 

The general search explores a six dimensional space: 
the location of the first circle centre, {9p, cf>p), the location of 
the second circle centre, {9r,())r), the angular radius of the 
circle a, and the relative phase of the two circles The 
number of operations needed to search for circle pairs scales 
as ~ A'^^logA''; thus implementation of a matched circles 
test is computationally very intensive. Because of this we 
will restrict our analysis to a search for pairs of circles cen- 
tered around antipodal points, so called back-to-back circles. 
Then, the search space can be reduced to four dimensions 
and the number of operations to ~ A'^'^ log A ^. We will not 
search for nearly back-to-back circles as in ICornish et al.l 
neither can we apply their hierarchical approach to 
speed up computations. In our case, we perform an analy- 
sis of high resolution maps (~ 20'), so that the approach of 
initially degrading them to a lower resolution and then refin- 
ing the search in higher resolution maps only for those circle 
pairs with the highest correlations risks missing some pairs 
of circles. As a consequence, our analysis is more computa- 
tional ly demanding than the analysis done bv I Cornish et al.1 
(|2004 )' preventing the study of a more general class of topolo- 
gies. Although the back-to-back search can detect a large 
class of the topologies that we might hope to find, there 
remain those topologies that predict matched circles with- 
out antipodal configurations such as the Hantzsche-Wendt 
space for flat universes or th e Picard space for hyperbolic 
universes (jAurich et al.ll2004l '). 

We used the healpix scheme with resolution parame- 
ter A'aide = 512 to define our search grid on the sky. To ac- 
commodate the use of the FFT approach, n — 2"^^^ points 
are used for each circle, where r corresponds to resolution 
parameter of the map given by A^sidc = 2^ . By definition 
this is also the angular resolution used to step through 0*. 
For a, we used steps a little bit smaller than characteristic 
scale 9c of coherent structures in the map, i.e. 2 The 
scale of coherence was approximated by the Full Width at 
Half Maximum (FWHM) of a gaussian beam profile used 
to smooth maps. The values of the temperature anisotropy 
at each point along the circle were interpolated based on 
values for the four nearest pixels with weights inversely pro- 
portional to the distances between a given point and the 
pixel centers. Other methods of interpolation can be used 
but we have verified that even using the temperature value 
of the nearest pixel does not change the value of the statistic 
significantly. 

To draw any conclusions from an analysis based on the 
statistic S'max(ct) it is very important to correctly estimate 
the threshold for a statistically significant match of the cir- 
cle pairs. The chance of random matches is inversely pro- 
portional to the number of coherent structures along the 
circles, therefore a false positive signal level of the statis- 
tic is especially large for circles with smaller radii. We used 
simulations of the maps with the same noise properties and 
smoothing scales as the WMAP data to establish the thresh- 
old such that fewer than 1 in 100 simulations would yield a 
false event. 
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It is important to note that the false detection level 
is smaller for higher resolution maps. Co nversely, as shown 
by ICornish. Spergel. fc StarkmanI ()l998l ) , the value of the 
statistic for matched circles 

depends on the ratio of the signal rms CTs to the noise rms 
ratio, denoted ^ for the map, thus the efficiency of the statis- 
tic to detect matched circles is increased by smoothing the 
data. A smoothing scale should therefore be adopted that is 
a reasonable trade-off between these two requirements. 

To eliminate the regions most contaminated by Galac- 
tic foreground residuals, we utilise the masks defined by the 
WMAP team. However, the statistic is very sensitive to the 
masking, particularly if a significant fraction of one or both 
circle pairs lies in the masked region. In this case, there is 
a significant probability to find a chance correlation of the 
temperature pattern between unmasked parts of the circles, 
thus increasing the false detection level. The effect is the 
most pronounced for circles with the largest radii, close to 
90 degrees, as well with very small radii. To avoid this, we 
restrict our analysis to those pairs for which less than half 
the length of each circle is masked. Though, the statistic 
computed in this way is not optimal, it is much more ro- 
bust with respect to masking. It should be noted that the 
dependence of the statistic on the cut could be avoided if 
the statistic were expressed as a function of the number of 
coherent structures along the circles. However, because our 
principal goal is to constrain the size of the Universe it is 
better to use the statistic determined as a function of circle 
radius. 

Finally, iKev et all l|200'iD have suggested that the per- 
formance of the statistic ([3| may be further improved by 
bandpass filtering the input map to minimise the anisotropy 
contribution from the ISW and Doppler terms which blur 
any signatures of topology present on the last scattering 
surface. We chose not to apply this kind of filtering because 
of possible complications arising from the interaction of the 
mask and the filter. Our analysis is enhanced by the use of 
high resolution maps rather than filtering. 



5 RESULTS 

Before beginning the search for pairs of matched circles in 
the WMAP data, we validated our algorithm using simula- 
tions of the CMB sky for a universe with 3-torus topology 
for which the dimension of the cubic fundamental domain 
L = 2c/ Ho, and with cosmolo gical parameters co rrespond- 
ing to the AC DPI model fsee lLarson et al.ll2010l . Table 3) 
determined from the 7- year WMAP results. In particular, we 
verified that our code is able to find all pairs of matched cir- 
cles in such map. The statistic S'^ax(o!) for the map is shown 
in Fig. [T] Note that the peak amplitudes in the statistic, 
corresponding to the temperature correlatio n for matched 
circles , decrease with radius of the circles. ICornish et al.l 
noted that this is primarily caused by the Doppler 
term which becomes increasingly anticorrelated for circles 
with radius smaller than 45°. The temperature match on 
the last scattering surface will be also diluted by the ISW 
effect which comes from the evolution of the structures close 



to the observer. However, to a large extent the m weighting 
used in the statistic (|3} mitigates against this and prevents 
the statistic from being dominated by the large scales. 

The intersection of the peaks in the matching statis- 
tic with the false detection level estimated for the W-band 
data (the same as in Fig. [Sjl defines the minimum radius of 
the correlated circles which can be detected for this map. 
The height of the peak with the smallest radius seen in Fig. 
[1] indicates that the minimum radius is about amin ~ 10°. 
To emphasise the advantage of using higher resolution data 
in the analysis, we also show the false detection level esti- 
mated for the ILC map (the same as in Fig. [2)). In this case, 
the minimum radius for detectable matched circles is about 

Qmin ~ 25°. 

The statistics for the WMAP ILC map are shown in 
Fig. [2] The map was analyz ed on both the fu ll sky and after 
applying the KQ85y7 mask l|Gold et al.ll2O10l ). In the case of 
the full sky analysis, the statistic shows some excess for 
pairs of circle with large radii. Nevertheless, after removing 
the most foreground contaminated regions the statistic does 
not reveal any unusually large values. This indicates that 
residuals of the Galactic emission, especially in the Galactic 
plane, cannot be neglected in the analysis of matched circles. 

It is interesting that similar behaviour is not seen for 
a full sky analysis of the first year WMAP ILC map. This 
probably explains why the importance of masking the most 
contaminated regions was not recognised in previous papers 
using the mat ched c i rcles statistic, the only exception being 
the paper by iThenI l|200d ). However, because they studied 
only the first year WMAP data, they also came to the con- 
clusion that the ILC map is good enough for the full sky 
analysis. We checked that the excess is not related to the bias 
correction (for residual foreground emission) applied to the 
7-year ILC map by the WMAP team since this was not im- 
plemented for the first year data. Indeed, the excess remains 
present for the 7-year ILC constructed by simply coadding 
th e individual f reque ncy channels with the weights provided 
in I Gold et al.l (I2OI0I ). Deeper studies of this problem are 
beyond the scope of this paper. Nevertheless, it provides a 
further warning against the naive use of a full sky ILC map 
for cosmological analysis, and particularly with respect to 
the two-point correlation function which is closely related 
to the statistic used her^ll. 

The false detection level shown in Fig. [2] was estimated 
on the basis of 100 simulations of the ILC map with the 
KQ85y7 mask applied, assuming a simple-connected uni- 
verse and established from the requirement that fewer than 
1 in 100 simulations should yield a false match. We did 
not correct the simulated maps for the bias coming from 
the random correlation between the CMB and the Galactic 
foregrouncjfl, since the effect is small outside of the Galactic 
plane region. Moreover, such a correction does not appear to 
be reliably evaluated as demonstrated with simulations that 
have been generated using the same foreground templates as 
for the simulated ILC maps. Maps corrected in this way do 



The remarks of lEfstathiou. Ma. fc Hansoij 1I2OIOI ') regarding the 
lack of evidence for Galactic foreground residuals at low Galactic 
latitudes are only relevant for the ILC map when smoothed to 
10° resolution. 

* so called 'Cosmic Covariance' 
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Figure 1. An example of the 5max statistic for a simulated CMB map of universe with the topology of a cubic 3-torus with dimensions 
L = 2 c/Hq. The dashed and dot-dashed lines show the false detection level for the statistic estimated from 100 MC simulations of the 
WMAP 7-year W-band map (the same as in Fig. 3) and ILC map (the same as in Fig. 2), respectively. 
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Figure 2. Smax statistic for the WMAP 7-year ILC map. The solid and dotted lines show the statistics Smax and Smax> respectively, 
for the ILC map masked with the KQ85y7 mask. The dash-dotted and dashed lines show the statistic Smax for the 7-year and the first 
year ILC unmasked maps, respectively. The false detection level estimated from 100 MC simulations of the ILC map and by scrambling 
the coefficients of the map are marked by dash-three dots and long dashed lines, respectively. The peak at 90° corresponds to a 
match between two copies of the same circle of radius 90° centered around two antipodal points. 
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not properly reflect uncertainties concerning details of the 
Galactic emission, especially in the Galactic plane. 

For comparison we show also the false detection level 
estimated on the basis of a 'scrambled' version of the ILC 
map (i.e. by randomly exchanging the m sphe rical harmonic 
coeffic ients, a^rn, of the map at every I) as in lCornish et al.l 
l|2004) . However, we still applied the KQ85y7 mask for each 
'scrambled' map to be consistent with simulations of the 
ILC data. As can be seen, the false detection level is slightly 
lower in this case. It is not surprising because the scrambling 
generates maps with the same two-point function and only 
different phase correlations that results in a smaller variance 
of the ^-statistic than for an ensemble of simulated ILC 
maps. 

In order to decrease the false detection level and be able 
to detect matched circles with smaller radius, we analyzed 
also the WMAP data with the highest angular resolution 
i.e. the W-band map, corrected for Galactic foregrounds and 
smoothed with a gaussian beam profile of FWHM = 20' to 
decrease the noise level. The statistic for this map analysed 
with the KQ85y7 mask is shown in Fig. [3] As for the ILC 
map, the false detection level was estimated from 100 simu- 
lations of the W-band data. 

We did not find any statistically significant correlation 
of circle pairs in either map. As seen in Fig. [T] the mini- 
mum radius at which the peaks expected for the matching 
statistic are larger than the false detection level is about 
Omin ~ 10° for the W-band map. Thus, we can exclude 
any topology that predicts matching pairs of back-to-back 
circles larger than this radius. This implies that in a fiat 
universe described otherwise by the best-fit 7-year WMAP 
cosmological parameters, a lower bound on the size of the 
fundamental domain is d = 2 i?Lss cos(ai„in) — 27.9 Gpc, 
where -Rlss is the distance to the last scattering surface. 

Of course, the above constraints do not apply to those 
universes for which the orientation of the matched circles 
is impossible to detect due to partial masking on the sky. 
Examples of topologies which can be overlooked are so 
called slab and chimney spaces (Riazuelo ct al. 2004a). For 
such topologies and appropriate configuration, all pairs of 
matched circles could lie in the Galactic plane that is re- 
moved by the mask. The probability of overlooking the cir- 
cles depends on the specific topology and radii of the cir- 
cles, so it is difficult to give a general expression. Neverthe- 
less, one can suppose that it decreases as the fraction of the 
masked sky decreases. 



6 CONCLUSIONS 

We have studied constraints on the topology of the Uni- 
verse using the method of matched circles as applied to the 
7-year WMAP ILC map and the foreground-reduced W- 
band map. We paid special attention to three aspects of the 
analysis that have been neglected in previous studies - the 
application of a mask, the use of high resolution data and 
the estimation of the false detection level on the basis of 
detailed MC simulations of the sky maps. 

The necessity for the application of a mask is due to the 
presence of residual Galactic foreground emission present 
even in the ILC map. These introduce a non-negligible effect 
on the matched circles statistic that is used for constraining 



topology. However, the possibility to apply the analysis to 
masked maps yields the opportunity to more tightly con- 
strain topology by using higher resolution, foreground cor- 
rected WMAP W-band data. Constraints on the topology 
depend significantly on the threshold for a significant match 
of the circle pairs. In order to estimate this correctly, we used 
100 MC simulations of the ILC and W-band maps assuming 
a simply-connected universe. The level of false detection cal- 
ibrated in this manner is slightly higher than that derived in 
ICornish et al] l|2004 ) using a method in which the analyzed 
maps are resampled by shuffling their spherical harmonic co- 
efficients. Although the difference is not very big, it is worth 
noting that the constraints on the size of Universe are over- 
estimated if one uses a lower level of false detection. 



The analysis of the WMAP W-band map, after correc- 
tion using templates of Galactic foreground emission, did 
not reveal any significant correlations for pairs of back-to- 
back circles with a radius greater than ~ 10° . This substan- 
tially extends the previous constraint on the minimum ra - 
dius of detectable matched circles given bv lKev et all (|2007| ) 
of 20° . It also places a lower bound on the size of fundamen- 
tal domain of about 27.9 Gpc for a fiat universe described 
by the best-fit 7-year WMAP cosmological parameters. Al- 
though this constraint concerns only those universes with 
such dimensions and orientation of the fundamental domain 
with respect to the mask that allow the detection of pairs of 
matched circles, the probability of overlooking circle pairs is 
rather low for the KQ85y7 mask that removes only a rela- 
tively small fraction of the sky. 



Of course, observations of the CMB with higher angular 
resolution and significantly lower noise level by the Planck 
satellite may yield even tighter constraints on the topology 
of the Universe. However, one should bear in mind that the 
possible improvement in the lower bound on the size of the 
Universe will not be substantial. The current constraint is 
not much smaller than the diameter of the observable Uni- 
verse 2 JiLSS = 28.3 Gpc, which imposes an upper limit on 
the size of the fundamental domain that it is possible to 
detect using the method of matched circles. The only sig- 
niflcant improvement might be related to improved mod- 
eling of the Galactic emission allowing the application of 
a smaller mask approaching closer to the Galactic plane. 
This would minimize the probability of overlooking topolo- 
gies with matched circles that are hidden within the masked 
region of the sky. 

Finally, as in ICornish et all (|2004l ). the studies could 
also be extended to search for nearly back-to-back circle 
pairs. However, the much higher computational require- 
ments for the analysis of high resolution maps make such 
studies difficult and extremely time-consuming. Neverthe- 
less, we can hope that the steadily increasing speed of pro- 
cessors and availability of larger computational resources 
will make such computations feasible in the coming years, 
thus allowing a flnal resolution of the problem of the topol- 
ogy of our Universe. 
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Figure 3. 5max statistic for the WMAP 7-year W-band map. Solid and dotted lines show the statistics Smax and 5max, respectively, 
for the W-band map masked with the KQ85y7 mask. The dashed line is the false detection level estimated from 100 MC simulations. 
The peak at 90° corresponds to a match between two copies of the same circle of radius 90° centered around two antipodal points. 
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